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t^- 1 ABSTRACT 



i—{ ■ Aims. We have studied the chemistry of nitrogen-bearing species during the initial stages of protostellar collapse, with a view 
t ' to explaining the observed longevity of N2H + and NH3 and the high levels of deuteration of these species. 

, Methods. We followed the chemical evolution of a medium comprising gas and dust as it underwent free-fall gravitational 
collapse. Chemical processes which determine the relative populations of the nuclear spin states of molecules and molecular ions 
were included explicitly, as were reactions which lead ultimately to the deuteration of the nitrogen-containing species N2H 4 " 
and NH3. The freeze-out of 'heavy' molecules on to dust grains was taken into account. 

Results. We found that the timescale required for the nitrogen-containing species to attain their steady-state values was much 
larger than the free-fall time and even comparable with the probable lifetime of the precursor molecular cloud. However, it 
transpires that the chemical evolution of the gas during gravitational collapse is insensitive to its initial composition. 
Conclusions. If we suppose that the grain-sticking probabilities of atomic nitrogen and atomic oxygen are both less than unity 
(S < 0.3), we find that the observed differential freeze-out of nitrogen- and carbon-bearing species can be reproduced by the 
model of free-fall collapse when a sufficiently large grain radius (a g ~ 0.50 (im) is adopted. Furthermore, the results of our 
collapse model are consistent with the high levels of deuteration of N2H + and NH3 which have been observed in L1544, for 
, example, providing that 0.5 < a g < 1.0 fxm. We note that the ortho:para H2D + ratio, and fractional abundance of ortho-H2D + , 
' which is the observed form of H2D + , should be largest where ND3 is most abundant. 
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1. Introduction 

Recent observational and theoretical studies of the initial 
stages of low-mass star formation have yielded some re- 
markable and intriguing results. The kinetic temperatures 
of the protostellar envelopes are observed to be very low 
and, as a consequence, atomic and molecular species con- 
taining elements heavier than helium ultimately freeze on 
to the grain surfaces (Bacmann et al. 2002, Tafalla et al. 
2004). As a result of this process, there remains a gas com- 
posed of atoms, molecules and molecular ions comprising 
H, D and He; this medium is, in effect, "completely de- 
pleted" of heavy elements. 

Observations show that essentially all C-containing 
species are depleted in the regions where CO appears to be 
absent (Tafalla et al. 2006); this is to be expected, as CO is 
the main gas-phase repository of carbon. The dissociative 
ionization of CO by He + produces C + , and reactions in- 
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volving C + yield C-containing molecules. Thus, it is very 
likely that the freeze-out of CO is accompanied by the de- 
pletion not only of directly related species, such as HCO + , 
but also of molecules such as HCN, H2CO and HC3N. On 
the other hand, N- or O-containing species such as OH, 
H2O, NH3 and N2H" 1 " are linked to the most abundant 
gas-phase forms of nitrogen and oxygen. 

In the case of nitrogen, it has long seemed likely that 
either N or N2 is the predominant gas-phase form; which 
of these species is more abundant is uncertain, in part be- 
cause the direct observation of either in dense clouds is not 
possible at the present time. Furthermore, chemical mod- 
els are bedevilled by large uncertainties in the rate coeffi- 
cients, at low temperatures, for several critical reactions, 
as will be mentioned below fSection l2.1J) . Nonetheless, it 
is clear that, in a situation where the freeze-out of species 
such as CO is occurring, the sticking probabilities (to 
grains) of N and N2 can become crucial for the nitrogen 
chemistry. 
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It is observed that the depletion of the heavy ele- 
ments on to the grains occurs differentially: the nitrogen- 
bearing species N 2 H + and NH 3 remain in the gas phase 
longer than molecules such as CO (Bergin & Langer 1997, 
Belloche & Andre 2004). In the context of a contract- 
ing cloud, "longer" implies higher densities. Tafalla et al. 
(2002) found that the fractional abundances of NH3 and 
N2H + in several prestellar cores either remain constant or 
increase with increasing density (see also Di Francesco et 
al. 2006). The relative longevity of NH 3 and N 2 H+ was 
believed to be due to a low binding energy of molecular 
nitrogen (their precursor) to the grains (Bergin & Langer 
1997, Aikawa et al. 2001). However, the recent laboratory 
work of Bisschop et al. (2006) has shown that the adsorp- 
tion energies of N 2 and CO are very similar. Furthermore, 
the probability that N 2 sticks when it collides with the 
surface of a grain at low temperature was measured and 
found to be close to unity. Thus, the longevity of N 2 H + is 
related neither to a low binding energy nor to a low stick- 
ing probability for N 2 , which may be expected to freeze 
out similarly to CO. In this case, the "volatile" nitrogen 
is most likely to be N. 

Another effect which is related to the depletion of 
heavy elements, owing to their adsorption by dust grain 
grains, is the enhancement of the abundances of the 
deuterated forms of those species remaining in the gas 
phase (see, for example, Ceccarelli et al. 2006). Given 
the relatively high abundances of N-containing species 
in prestellar cores, it is perhaps not surprising that the 
deuterated forms of ammonia are much more abundant 
than anticipated on the basis of the elemental ud / nu ra- 
tio (cf. Roueff et al. 2005). In the present paper, we in- 
terpret these observations in the contexti of a model of a 
core collapsing on the free-fall timescale. In order to re- 
produce the observations, it is necessary that the timescale 
for deuterium fractionation should be comparable to both 
the dynamical timescale and the timescale for depletion on 
to dust grains; this requires (n g /n-ft)ira,g, where n g is the 
number density and a g is the radius of the grains, to be an 
order of magnitude smaller than is observed in the diffuse 
interstellar medium (Mathis et al. 1977), or, equivalently, 
the mean radius of the grains to be an order of magnitude 
larger. 

In Section 2 we discuss the relevant chemistry and, in 
an Appendix, the treatment of the ortho and para forms 
of species such as NH3 . Section 3 contains our results and 
their comparison with the relevant observations, particu- 
larly of the prestellar core L1544. In Section 4, we make 
our concluding remarks. 



2. Chemistry 

2.1. Nitrogen-containing species 

The nitrogen-containing species which are observed in 
prestellar cores are N 2 H + and NH3. However, elemental 
nitrogen is expected to be mainly in the form of N and 



N 2 , in the gas phase. The principal reactions converting 
N to N 2 are 

N + OH^NO + H (1) 
and 

N + NO -> N 2 + O (2) 

(Pineau des Forets et al. 1990). Subsequent reactions of 
N 2 (i) with H+ yield N 2 H+ and (ii) with Hc+ yield N+, 
which initiates a series of hydrogen-abstraction reactions 
with H 2 , ending with NH^, which dissociatively recom- 
bines with electrons, yielding NH3. 

Pineau des Forets et al. (1990) assumed that the re- 
actions (JTJ and O of N with OH and NO with N had 
small (50 K) barriers. At the kinetic temperature consid- 
ered here (T = 10 K), the barriers contribute to making 
these reactions very slow: using the equilibrium fractional 
abundances of OH or NO, the timescale characterizing the 
conversion of N to N 2 is in excess of 10 s yr in molecular 
gas of density tt-h — 10 4 cm" 3 . Under these circumstances, 
the assumption of static equilibrium, when determining 
the initial composition of the gas, is unlikely to be valid 
in so far as N 2 and the other products of the nitrogen 
chemistry are concerned. We shall return to this point in 
Section 13.11 

To date, there has been no independent experimental 
or theoretical verification of the existence of barriers to the 
reactions Q and J5J) of N with OH and NO. Furthermore, 
as these are 'radical-radical' reactions (i.e. reactions be- 
tween two open-shell species), it is possible that there are 
no such barriers (Smith 1988). Accordingly, in the present 
study, we have removed the barriers to these reactions; 
at T = 10 K, this choice has the consequence of enhanc- 
ing their rate coefficients, and reducing the corresponding 
timescales, by a factor of 150, to the order of 10 6 yr; the as- 
sumption that the initial abundances are those computed 
in steady state then becomes more reasonable. However, 
these reactions remain slow in the context of the free- 
fall collapse, considered below, for which the timescale is 
4.3 x 10 5 yr when — 10 4 cm -3 initially. It follows that 
simulations of the collapse of the core by means of calcu- 
lations which assume that equilibrium is attained at each 
(increasing) value of njj will not predict the evolution of 
the nitrogen chemistry correctly. 

Molecular nitrogen (N 2 ) is destroyed in the reactions 



N 2 + H+ -> N 2 H+ + H 2 (3) 

N 2 + Hc+ -> N+ + He (4) 
and 

N 2 + Hc+ -> N+ + N + He (5) 



The rate coefficients for these reactions are compara- 
ble and are of the order of the Langevin rate coefficient 
(10 -9 cm 3 s _1 ). In static equilibrium at nn = 10 4 cm -3 , 
the fractional abundance of H3 , n(Hg )/nn ~ 2 x 10~ 9 , 
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is an order of magnitude larger than the fractional abun- 
dance of He + , 7i(He + )/n H ~ 2xl0~ 10 . Thus, the timescale 
for destroying N 2 is initially of the order of 10 6 yr, com- 
parable with the timescale for its formation and also with 
the timescale for its freeze-out. The total fractional abun- 
dance of H;j~ and its deuterated forms and the fractional 
abundance of He + remain approximately constant during 
the free-fall collapse, up to a density tt-h = 10 7 cm" 3 . 
On the other hand, the fractional abundance of OH falls 
rapidly, owing to freeze-out, for n-n ;> 10 5 cm" 3 . It fol- 
lows that N2 is removed not only by freeze-out but also 
by gas-phase chemical reactions during the early phases 
of collapse, when it ceases to be formed. As N2 is the 
precursor of both N2H + and NH3, the abundances of the 
latter species during collapse are dependent on the initial 
abundance of N2 . 

The main process of destruction of N 2 H + is dissocia- 
tive recombination: 

N 2 H+ + c- -> N 2 + H (6) 

N 2 H+ + e" -> NH + N (7) 

Recent ion-storage-ring measurements (Geppert et al. 
2004) have shown that, somewhat counter-intuitively, the 
branching ratio to the products of reaction © is 0.36, 
smaller than the branching ratio of 0.64 to the products 
of reaction @ • We have adopted the experimental results 
of Geppert et al., including their expression for the total 
dissociative recombination coefficient [summed over the 
products of © and ©], namely 1.0 x 10" 7 (T c /300)-°- 51 
cm s _1 , where T c is the electron temperature. 

2.2. Oxygen-containing species 

Specific mention should be made also of some of the un- 
certainties in the oxygen chemistry. Neither H 2 nor OH 
is readily observable in cold dense cores: H 2 lacks easily 
excitable emission lines, and the the 18 cm lines of OH can 
presently be observed only at angular resolutions such that 
any emission from the core is swamped by that from a less 
dense but warmer envelope. Consequently, we have no di- 
rect observational evidence of the survival of species such 
as OH, H 2 and H 3 + in the densest regions (nn > 10 5 
cm" 3 ) of prestellar cores. However, in somewhat less dense 
regions, SWAS data indicate that neither water nor molec- 
ular oxygen is a major repository of oxygen in the gas 
phase (Bergin et al. 2000; see also Lee et al. 2004). These 
authors conclude also that the best agreement with the 
SWAS observations is obtained if the gas-phase C:0 ele- 
mental abundance ratio is greater than 0.9, in which case 
CO contains most of the gas-phase oxygen. This situa- 
tion could arise if oxygen is mainly in the form of water 
ice at densities for which CO is still in the gas phase. We 
have borne this result in mind when specifying the initial 
conditions of some of our models (Section 
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2.3. Deuteration reactions 

Under conditions of complete heavy element depletion, the 
principal molecule is H2, and He remains in atomic form. 
The main ions are H + and H^, together with its deuter- 
ated forms H 2 D+, D 2 H+ and (Roberts et al. 2003, 
2004; Roueff et al. 2005). Both H 2 D+ and D 2 H+ have been 
observed in prestellar cores (Caselli et al. 2003, Vastel et 
al. 2004); is more difficult to observe, owing to the 
absence of a permanent dipole moment. The remarkably 
high levels of deuteration of H^" which are observed arise 
from the differences in the zero-point vibrational energies 
of ions in the deuteration sequence, for example, between 
H^" and H2D + , which has a larger reduced mass. At low 
kinetic temperatures, the corresponding deuteration reac- 
tion 

H+ + HD ^ H 2 D+ + H 2 (8) 

is strongly favoured, energetically, in the forwards direc- 
tion, when ground state reactants and products are in- 
volved. The other reactant, HD, is formed, like H 2 , on the 
surfaces of dust grains. 

In earlier papers (Walmsley et al. 2004, Flower et al. 
2004, 2005, 2006), we considered the deuteration of H| 
and emphasized the importance of distinguishing between 
the ortho and para forms of the reactants and the prod- 
ucts. Different nuclear spin states are, through the appro- 
priate quantum statistics (Fermi-Dirac or Bose-Einstein), 
associated with different rotational states of the molecule. 
In the best-known example of H 2 , ortho states (with total 
nuclear spin I — 1) are associated with rotational levels 
with odd values of J, whereas para states (I = 0) are 
associated with even values of J. Thus, the lowest ortho 
state, J = 1, lies approximately 170 K higher than the 
lowest para state, J = 0, which is the ground state of 
the molecule. This energy, which is large compared with 
the thermal energy of the gas (the kinetic temperature 
T « 10 K) becomes available in the backwards reac- 
tion jHJ, in the case where ortho-H2 is the reactant with 
H2D+. As the backwards reaction reverses the deuteration 
of , it is essential to know the ortho:para H2 ratio in or- 
der to predict correctly the degree of deuteration (Gerlich 
et al. 2002, Flower et al. 2006). 

It is well established that H^ plays a pivotal role 
in interstellar chemistry. Similarly, its deuterated forms, 
H2D+, D 2 H + and , are determinants of the degree of 
deuteration of other species. The reaction 

N 2 + H 2 D+ # N 2 D+ + H 2 (9) 

and the analogous reactions with D2H + and Djj~ are the 
keys to understanding the N 2 D + :N 2 H + ratio in prestel- 
lar cores. The degree of deuteration of NH3 is determined 
by analogous reactions with the deuterated forms of H^ . 
Thus, the interpretation of the (high) levels of deutera- 
tion of both N 2 H + and NH3 observed in prestellar cores 
involves consideration of the ortho:para abundance ratios 
of certain gas-phase species. 
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3. Results 

We have considered the gravitational contraction (free- 
fall collapse) of a condensation of gas and dust with an 
initial number density, riu = n(H) + 2n(H2), in the range 
10 3 < tih < 10 5 cm~ 3 and kinetic temperature T — 10 K. 
The initial value of the grain radius was varied in the 
range 0.05 < a g < 1.0 /im, for reasons that will become 
apparent below. The smallest value, a g = 0.05 /im, yields 
approximately the same grain opacity, n g Tra^, as the grain 
size distribution of Mathis et al. (1977), with limits of 
0.01 < a g < 0.3 fj,m, and where n s is the grain number 
density. 

The initial composition of the gas was computed on 
the assumption of static equilibrium ("steady state"). 
The elemental abundances and initial depletions are 
specified in table 1 of Flower et al. (2005); they yield 
an initial gas-phase C:0 elemental abundance ratio, 
nc/no = 0.67, and a dust:gas mass ratio of 0.0094. 
In the light of the discussion in Section 12.21 higher 
values of uq/uo were also considered. The chemistry 
comprised reactions involving species containing H, He, 
C, N, O and S and distinguished between the possible 
nuclear spin states of and its deuterated forms, 
and those of H^ and H2. It also distinguished be- 
tween the nuclear spin states of nitrogen-containing 
species, following the discussion in Appendix El 
The chemistry and list of species are available from 
http:// massey.dur.ac.uk/drf/ protostellar /chemistry .species 

In order to study the deuteration of nitrogen- 
containing molecules, such as NH 3 , and ions, such as 
N2H + , we included the relevant chemical reactions but did 
not introduce explicitly the different nuclear spin states 
of the deuterated forms. Thus, NH2D, for example, was 
treated as a single chemical species; we did not distinguish 
between its ortho and para forms. This simplification was 
necessary to ensure that the chemistry remained manage- 
able, but it did not, in our view, compromise the reliabil- 
ity of the results relating to the degrees of deuteration of 
nitrogen-containing species. 

In most of the models, the grain sticking probability 
S(N) — 0.1 for atomic nitrogen, and S(0) — 0.1 also 
(see the discussion in Section l3~2|> . Larger values of 5*(N) 
and S(0) were considered, and reference to the results 
of these calculations will be made below. For all other 
atomic and molecular species, S — 1.0. The rate of cosmic 
ray ionization of hydrogen was taken to be C = 1 X 10~ 17 
s _1 , as in our earlier papers (Flower et al. 2005, 2006). 

3.1. Initial abundances 

As we noted in Section the timescales associated with 
the chemistry of the nitrogen-bearing species are partic- 
ularly large - even when, as here, we assume that there 
are no barriers to the initiating reactions and ©• To 
illustrate this point, we plot, in Fig.^ the fractional abun- 
dances of gas-phase N, N2, NO and NH3 as functions of 
time, t, at a constant gas density of Uh = 10 4 cm -3 , start- 



ing with the nitrogen in atomic form. Freeze-out of the 
heavy elements is neglected in this calculation. 

Figure ^ shows that the formation of molecular nitro- 
gen requires t«5x 10 6 yr; the fractional abundances of 
NO and NH3 attain steady state on a similar timescale, 
as they are a precursor to and a consequence of, respec- 
tively, the formation of N2 . This timescale is comparable 
with the lifetimes, r, of molecular clouds, as derived by 
Tassis & Mouschovias (2004; r w 10 7 yr), but greater than 
the ages of clouds in the solar vicinity, as determined by 
Hartmann et al. (2001; r ~ 10 6 yr). Thus, it seems pos- 
sible, though by no means certain, that such species have 
reached their steady-state abundances before protostellar 
collapse begins. 

In fact, the nitrogen chemistry is peculiar, in that it 
is activated by a reaction between two neutral species 
(N and OH) of low fractional abundance (of the order 
of 10 -5 and 10 -8 , respectively, in steady-state). Overall, 
the steady-state chemistry is dominated by ion-neutral 
reactions involving H2 as one of the reactants. Although 
the timescale for cosmic ray ionization is very long, of the 
order of w 10 9 yr, the principal ionization processes 
involve H2 and He, which are the most abundant species 
in the gas; the charge is then passed to 'heavy' species, 
in rapid ion-neutral reactions. The resulting rate of ion- 
ization of the heavy species, X, is approximately nn/^x 
faster than the rate of direct cosmic ray ionization of X. 
As a consequence, the abundances of ions and the species 
which form through ion-neutral reactions, such as OH, 
reach steady state on a timescale which is about an order 
of magnitude smaller than the free-fall time of 4.3 x 10 5 yr 
for an initial density «h = 10 4 cm~ 3 . 

In view of the uncertainty in the gas-phase C:0 ele- 
mental abundance ratio (see Section l2~2"|) . we recalculated 
the steady-state composition of the gas with a higher 
value of this ratio, no /no = 0.95. We found that, although 
the fractional abundances of OH and NO decreased, by a 
factor of aproximately 2.5, that of N2 fell by only 15%, 
with similar decreases for N2H + and NH3. We conclude 
that, as long as the gas-phase C:0 abundance ratio re- 
mains less than 1, there remains sufficient oxygen in the 
gas phase in the form of OH to drive the nitrogen chem- 
istry. In any case, reactions such as 



He + 
and 



CO -> C+ + O + He 



He+ + 2 -> 0+ + O + He 



(10) 



(11) 



restore atomic oxygen to the gas phase, thereby enabling 
OH to form in the sequence of reactions initiated by H^ (O, 
H 2 )OH+. 

We note that the analogous reaction JSJ of He + with 
N2 not only leads ultimately to the formation of NH3 but 
also restores atomic nitrogen to the gas phase. For this 
reason, the chemical evolution predicted by the model cal- 
culations proves to be insensitive to the initial value of the 
atomic to molecular nitrogen abundance ratio; see the re- 
sults in the following Section I3~2l 
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3.2.1. Dependence on the initial conditions 
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Fig. 1. The fractional abundances, relative to ~ 
n(H) + 2n(H 2 ), of N, N 2 , NO and NH 3 as functions of 
time, t, through to the attainment of their steady-state 
values. A constant gas density of njj = 10 4 cm -3 and 
temperature T — 10 K were assumed and freeze-out of 
the heavy elements was neglected. 

3.2. Free-fall collapse 



The free-fall time is given by 



3tt 



32G/9Q 



where po is the initial mass density; Tg = 4.3 x 10 5 yr for 
an initial number density nn = 10 4 cm~ 3 . We have as- 
sumed that the density evolves on the free-fall timescale. 
(We note that Ward-Thompson et al. (2006) have pre- 
sented recently statistical arguments which suggest that 
core lifetimes are a few times Tff, possibly owing to mag- 
netic support). 

The freeze-out time 



Tfo = [ngTTdgVthS] 



where n g 7rOg is the grain opacity and Vth is the thermal 
speed of the species which sticks to the grain, with a stick- 
ing probability S. For a g = 0.05 pm, a molecule such as 
CO freezes out on a timescale Tf Q = 2.2 x 10 5 yr when 
riH = 10 4 cm~ 3 . Thus, the free-fall and the freeze-out 
times are comparable at this density and for this value 



of the grain radius. We note that n g 



,-3 



and hence 



Tf oc a g : freeze-out occurs at later times and consequently 
higher densities if the grains are larger; this effect will be 
seen in the results presented below. 



The uncertainties in the initial chemical composition of 
the gas has already been considered on Section |3~T1 and, 
in particular, the fact that the chemistry may not have 
attained steady state before gravitational collapse begins. 
Accordingly, we have carried out a number of test calcu- 
lations, varying the intial composition from equilibrium. 
These calculations were designed to test the sensitivity of 
the fractional abundance profiles, during the subsequent 
free-fall collapse, to such variations. In Fig. |2 is shown 
a sample of such calculations. Fig. [5] enables the conse- 
quences of assuming that nitrogen is initially in atomic 
rather than mainly in molecular form to be appreciated. 
There is a large effect initially on the NH3 and N2H" 1 " 
abundances, but, by densities of a few times 10 5 cm~ 3 , 
the differences have largely disappeared. Furthermore, and 
for reasons given above, the effect of setting nc/no = 1-0 
initially, and thus reducing drastically the oxygen avail- 
able to form OH, is less than might have been antici- 
pated, although the formation of N2 is slower than when 
nc/no = 0.67. 

Fig-Eldemonstrates the remarkable insensitivity of the 
chemical evolution during collapse to the initial condi- 
tions, even in the case of the nitrogen-containing species, 
whose chemistry involves relatively long timescales (cf. 
Section 12. Major variations of the initial conditions 
from steady-state, such as assuming that all the gas-phase 
nitrogen is initially in atomic form, have transient and 
limited effects on the subsequent evolution. In essence, 
Fig. |21 confirms the effectiveness of the ion-neutral chem- 
istry, which we have already considered in Section l3.ll 

3.2.2. Differential depletion on to grains 

As noted in the Introduction, an intriguing observational 
result, which awaits an interpretation, is the relative 
longevity of N2H 4 " and NH3 during freeze-out of the heavy 
elements from the gas phase. Explanations in terms of a 
low adsorption energy or a low sticking probability for 
N2, the precursor of both N 2 H + and NH3, appear to have 
been invalidated by recent experimental work (Bisschop et 
al. 2006). However, we are unaware of any measurements 
of the sticking probability, S'(N), of atomic nitrogen (or 
S(0), which is also relevant to the gas-phase chemistry; 
see below). N and O are lighter than N2 and have corre- 
spondingly higher vibrational frequencies when interact- 
ing with the grain surface, which, in the present context, 
is probably composed essentially of CO, deposited during 
freeze-out. The issue which waits to be addressed, exper- 
imentally or theoretically, is the strength of CO - N and 
CO - O bonds, relative to CO - CO (or CO - N 2 ). We 
note that there exists some experimental evidence for the 
enrichment of CO, relative to N 2 , during freeze-out at low 
temperatures, although on to a layer of predominantly wa- 
ter ice, rather than CO ice (Notesco & Bar-Nun 1996). 

The formation of molecular nitrogen is initiated by 
reaction Q), N(OH, H)NO, which requires both N and 
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OH to be available in the gas phase. The hydroxyl rad- 
ical, OH, is formed in the dissociative recombinations of 
H 2 + and H 3 + , which are produced from atomic oxy- 
gen, through an initial proton-transfer reaction with H^ 
and sub sequent hydrogenation of OH + and H20 + by H2. 
These ! ;hemical considerations suggest that relatively low 
sticxin 1 probabilities for both atomic nitrogen and atomic 
oxyger might contribute to the longevity of N2H 4 " and 
NH3. iv low value of 5(0) would also have consequences 
for thd abundances of 'pure' oxygen-containing species, 
such a:s O, O2 and H 2 0; but the angular resolution and 
sensitivity of current instrumentation are such that the 
core cannot be distinguished from the more extended, 
lower density regions that are observed, for example, in 
OH (Goldsmith & Li 2005). 

Although our results are not directly comparable to 
u\h observations (which have a linear spatial resolution 
"vr4Njly in the range 0.15-0.3 pc), it is unlikely that the 
cn^rnS^al abundances at densities of order 10 4 cm -3 (our 



Fig. 2. The fractional abundances, relative to nn ~ 
n(H) + 2n(H 2 ), of (a) NH 3 and (b) N 2 H+, as functions 
of riH, during free-fall gravitational collapse from an ini- 
tial density nn = 10 4 cm -3 . A constant temperature 
T = 10 K and a constant grain radius a g = 0.20 /itm 
were assumed. The full curves correspond to an elemental 
C:0 ratio, in the gas phase, of nc/no = 0.67 (our stan- 
dard assumption), and the broken curves to nc/no =1.0. 
In the "molecular" case, the gas was initially in molec- 
ular form (our standard assumption), whereas, in the 



initialj ^yaluc) differ greatly from those observed in clouds 
with densities of a few times 10 3 cm -3 , which are believed 
to be relevant to the measurements of OH, for example. 
Goldsmith & Li (2005) found n(OH)/n(H 2 ) = 2 x 10~ 8 
in L1544, Bergin & Snell (2002) obtained upper limits to 
n(ortho-H 2 0/n(H 2 ) of the order of 3 x 10~ 8 in two clouds, 
and Pagani et al. (2003) derived n(0 2 )/n(H 2 ) < 10~ 7 
in several sources. The observational results for O2, in 




particular, plead in favour of models in which nc/no = 
1.6 mi lally in the gas phase: these models have steady- 
state C 2 abundances close to the observed upper limit. On 
the other hand, when nc/no — 0.67 initially, the steady- 
state abundance of O2 is roughly two orders of magnitude 
larger. However, we recall that the nitrogen chemistry is 
not sensitive to the initial value of nc/no- 

In Fig.|3arc shown the profiles of N2, NH3 and N2H + , 
as functions of nn, in the course of free-fall collapse from 
an initial steady-state composition and an initial density 
fiH =-10 4 cm -3 ; the profile of CO is also plotted, for 
\\cchtrpa' ison. The values adopted for 5(N) = 5(0) are in- 
atcY in panel (a) of Fig. |31 and 5 = 1.0 for all other 
Qtomk; and molecular species. As 5(N) and 5(0) are re- 
"Succc^Qatomic nitrogen and also OH, which is produced 
from atomic oxygen, are removed from the gas phase more 
slowly than CO, and the nitrogen chemistry continues 
to be driven by reactions JU and (J2J. Consequently, the 
abundances of the nitrogen-containing species increase, 
relative to CO, as the density approaches nn = 10 6 cm -3 . 
Supplementary calculations have confirmed that the in- 
crease in the N2H + :CO ratio, for example, occurs only 
when 5(N) and 5(0) are reduced simultaneously. Thus, 
we conclude tentatively that the longevity of N2H + and 
NH3 is related to low grain sticking probabilities for both 
atomic nitrogen and oxygen. 

We show the effect of varying the grain size in Fig. 
The onset of freeze-out is delayed to higher gas densi- 
ties as a g increases, owing to the reduction in the to- 
tal grain surface area, which, in our model, is inversely 
proportional to a g . Thus, a g = 0.5 /im corresponds to 
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a value of (n g /nn)TTa^ = 1.7 x 10 cm , which is an 
order of magnitude smaller than is deduced from obser- 
vations of dust in the diffuse interstellar medium (see, for 
example, Mathis et al. 1977). CO is depleted by a factor 
of e by a density njj ~ 6 x 10 4 cm~ 3 for a g = 0.20 /j,m 
and riH ~ 2 x 10 5 cm~ 3 for a g = 0.50 /urn. These val- 
ues of nH may be compared with the densities char- 
acterizing CO freeze-out, deduced from observations of 
five nearby cores (Tafalla et al. 2002, table 4), namely 
3 x 10 4 <; n H <; 1 x 10 5 cm" 3 (note that Tafalla et al. 
expressed their results in terms of n(H2) « tih/2). 

It is interesting that, according to our models (Fig.[3J|, 
there is a wide range of density, 1 x 10 5 < rtn <^ 2 x 
10 6 cm" 3 , in which the fractional abundance of N2H + 
is roughly constant. After an initial drop, the fractional 
abundance of ammonia goes through a maximum at nn ~ 
10 6 cm -3 ; this is reminiscent of the deduction of Tafalla et 
al. (2002), that, in the five cores which they observed, the 
fractional abundance of N2H + remains roughly density- 
independent whilst that of NH3 appears to increase with 
increasing core density. Taking into account the uncertain- 
ties in the chemistry, we find this qualitative agreement to 
be encouraging. 

With the comments in the previous paragraphs in 
mind, we display, in the upper panels of Fig. [S] the frac- 
tional abundances of CO, N 2 H+ and NH 3 , observed in 
the prestellar core L1544 (Tafalla et al. 2004), along with 
the estimated error bars; the density, nn, has been de- 
duced from the peak intensity of the dust emission. The 
abundance profiles which derive from the free-fall model, 
with an initial density nu = 10 4 cm" 3 and grain radii 
a g = 0.5 /im and a g = 1.0 /im, are also plotted in Fig. El 
The observed fractional abundances are reproduced sat- 
isfactorily by the model with a g = 1.0 /im. As we shall 
see in Section EOl below, the evidence from the calculated 
and observed levels of deuterium fractionation of N2H + 
and NH3 is also consistent with this value of a g . 



3.3. Deuteration of nitrogen-containing species 

Observations of prestellar cores have shown that, not only 
are nitrogen-containing species relatively long lived, but 
also they exhibit remarkably high degrees of deuteration. 
Indeed, even triply-deuterated ammonia, ND3, has been 
detected (van der Tak et al. 2002, Lis et al. 2002), with an 
abundance, relative to NH3, which is more than 10 orders 
of magnitude larger than would be anticipated on the basis 
of the tid/^h (~ 10 -5 ) elemental abundance ratio. 

We have already considered, in Sect ion l2~3l the mecha- 
nism of deuteration of nitrogen-containing species such as 
NH3 and N2H+, which involves reactions with the deuter- 
ated forms of Hj. H2D+ is formed through the deutera- 
tion of H3" by HD and (together with D2H + and ) is 
the main agent of deuteration of other species. The frac- 
tional abundances of the deuterated forms of become 
so large that the timescale characterizing the deuteration 
of N 2 H + , in reaction (JSJ and the analogous reactions with 
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Fig. 3. The fractional abundances, relative to nn, of (a) 
N2, and (b) CO, NH3 and N2H" 1 ", as functions of nn, 
predicted by the free-fall model with an initial density 
nn = 10 4 cm -3 and for a grain radius a g = 0.5 /im. The 
adopted values of the sticking probabilities 5(N) = 5(0) 
are indicated in panel (a), and the key is identical in panel 
(b). 
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D 2 H + and D3 , falls to values comparable with the free- 
fall time. Consequently, the level of deuteration of N 2 H + 
increases with the gas density, in the course of the 
collapse. 

As shown by Roueff et al. (2005, Table 8; see also Lis et 
al. 2006), the population densities of successive stages of 
deuteration of NH3 are observed to decrease by only about 
an order of magnitude, and 10 -4 < n(ND3)/n(NH3) < 
10 -2 in the prestellar cores which have been observed. 
Similarly, N 2 D + is observed to be approximately an or- 
der of magnitude less abundant than N 2 H + (Caselli et 
al. 2002). In Fig. these isotopic abundance ratios are 
plotted, together with the fractional abundances of NH 3 
and its deuterated forms. The fractional abundances of 
the deuterated forms of ammonia attain a maximum in a 
rather narrow range of density 5 x 10 5 < nn 5 x 10 6 
cm~ 3 . The maximum fractional abundance shifts to higher 
densities with increasing degree of deuteration and with 
increasing grain size. Thus, ND3, for example, should sam- 
ple higher density gas than NH 3 . 

Also plotted in Fig. is the ortho:para H 2 ratio, in 
view of the key role which it plays in the deuteration pro- 
cess. The ortho:para H 2 ratio tends towards thermaliza- 
tion later (i.e. at higher rtu) with increasing grain size. 
Because freeze-out occurs later for larger a g , the fractional 
abundances of H^" and H + maximize later, which, in turn, 
delays ortho-to-para conversion in proton exchange reac- 
tions of these species with H 2 ; this has knock-on effects 
on the ortho:para ratios of other species, such as H 2 D + . 

The behaviour of H^" and its deuterated forms is shown 
in Fig. [7| for a g = 0.5 fim. The increase in the degree of 
deuteration of , as freeze-out occurs, is apparent in 
this Figure; H 2 D + , D 2 H + and are responsible for the 
deuteration of NH 3 . The rapid decrease in the abundance 
ratios of the ortho and para forms of H^" and H 2 D + is 
seen, in Fig. [7jp, to coincide with the decrease in the or- 
tho:para H 2 ratio; these ortho:para ratios are linked chem- 
ically, through proton-exchange reactions (cf. Flower et 
al. 2006). The ortho:para H 2 D + ratio depends on both 
nH and a g . The maximum fractional abundance of ortho- 
H 2 D + is attained for rtn ~ 10 6 cm -3 and is of the order 
of 10~ 10 . The maximum in ortho-H 2 D + occurs close to 
the maxima in the deuterated forms of NH3, suggesting 
that observations of these species sample the same region 
of the core. The fact that the ortho:para H 2 D + ratio re- 
mains large for densities up to nn ~ 10 6 cm -3 probably 
accounts, at least partly, for the ease with which ortho- 
H 2 D + is being detected in prestellar cores. 

In Fig.0 we compare the predictions of our model with 
the observations of the degrees of deuteration of N 2 H + and 
NH 3 in the prestellar core L1544. The observational and 
theoretical results are in satisfactory agreement for cither 
of the values of a g for which calculations are displayed, 
i.e. a g = 0.5 /im or a g = 1.0 /iin. Grain sizes a g < 0.5 /^m 
or a g > 1.0 /urn would appear to be excluded by the ob- 
servations, when interpreted in the context of our model. 



3.4. Relative population densities of nuclear spin states 

We note first that our expectations regarding the rela- 
tive populations of the states I = 0, J = 1 and 1 = 2 
of NH4 (Appendix lA.ljl and the para:ortho NH3 ratio 
(Appendix IA.2|I are confirmed by the results of the free- 
fall collapse model, which show these population ratios 
to be independent of the ortho:para H 2 ratio and of the 
gas density and equal to the numerical values derived in 
Appendix lAl 

As noted in Appendix I A. 31 the only species for which 
more than one nuclear spin state has been observed in 
prestellar cores is NH 2 D. The ortho:para NH 2 D ratio 
which is predicted, on the basis of the arguments in 
Appendix lA.3l is 2:1; this may be compared with the value 
of 3:1 which is anticipated from the ratio of the degenera- 
cies, 21+ 1, of the corresponding nuclear spin states, 1 = 1 
and / = 0. In practice, the observations reported by Shah 
& Wootten (2001) are perhaps marginally more consistent 
with a ratio of 3:1 than 2:1; but, allowing for the associ- 
ated observational error bars, the lower value (of 2:1) is 
consistent with their observations of 4 out of a total of 8 
sources. 

4. Concluding remarks 

— We have investigated the chemistry of nitrogen- 
bearing species during the initial phases of protostel- 
lar collapse. We find that the timescales for these 
species to reach their steady-state, equilibrium values 
are large - comparable to, or perhaps greater than, the 
lifetimes of the precursor molecular clouds. Even if, as 
we have supposed, steady state is attained prior to the 
onset of gravitational collapse, the nitrogen chemistry 
departs from steady state as soon as collapse begins. 
Large discrepancies develop, as rh increases in the 
range 10 4 < nn < 10 7 cm -3 , between the predictions 
of the free-fall model and the corresponding results in 
steady-state, for important nitrogen-bearing species, 
specifically N 2 , NO, NH 3 and N 2 H+. The steady-state 
calculations underestimate the fractional abundances 
of N 2 and N 2 H + , and overestimate those of NO and 
NH 3 , by factors which approach 2 to 3 orders of mag- 
nitude in the range of gas densities considered here. 

— We have assumed that the collapse occurs rapidly, on 
a timescale of approximately 4 x 10 5 years. A longer 
evolutionary timescale would allow CO to freeze-out 
at lower densities than are observed. In addition, this 
timescale is compatible with current observational es- 
timates of the lifetimes of prestellar cores. Another as- 
sumption of our model is that the initial chemical com- 
position corresponds to the gas phase being in steady 
state at the initial density. In fact, we have found that 
the chemical evolution predicted by our model is in- 
sensitive to this assumption, and, in particular, to the 
n c l n o elemental abundance ratio and even to the ex- 
tent to which elements heavier than helium have been 
incorporated into molecules and molecular ions. We 
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have made ad hoc assumptions about the sticking co- 
efficients, S, of atomic nitrogen and atomic oxygen. 
Clearly, laboratory work is needed to check the values 
of S(0) and 5(N) which have been adopted. 

— We have considered the issue of the relative longevity 
of nitrogen-containing species: N 2 H + and NH3 are 
observed to be still present in the gas phase when 
other molecules, specifically carbon-containing species 
such as CO, have already frozen on to the grains. 
These observations suggest that a major reservoir of 
elemental nitrogen, probably N or N2, has either a 
low adsorption energy or a low probability of stick- 
ing to grains. Recent experimental work (Bisschop et 
al. 2006) excludes these possibilities in the case of N 2 . 
Accordingly, we have supposed that the sticking prob- 
ability for atomic nitrogen, S(N) < 1. Because the 
chemical pathways leading to N 2 H + and NH3 are ini- 
tiated by N(OH, H)NO [reaction (jTJ above], and OH 
is produced from atomic oxygen, chemical consider- 
ations imply that 5(0) < 1 also. Subject to these 
hypotheses, we find that the differential freeze-out of 
nitrogen- and carbon-bearing species can be repro- 
duced by our model of free fall collapse when a suf- 
ficiently large grain radius (a g « 0.5 /im) is adopted. 
Although OH cannot be observed with sufficient angu- 
lar resolution to investigate directly its chemical role 
in prestellar cores, NO can and has been observed, in 
L1544 and L183; these observations will be presented 
in a forthcoming publication (Akyilmaz et al. 2006, in 
preparation) . 

— In effect, we have varied the values of the timescales 
characterizing the gas phase chemistry and depletion 
on to grain surfaces, relative to the timescale for free- 
fall collapse, investigating the consequences for observ- 
able quantities, such as the degree of deuteration of 
ammonia. Our results indicate that (n g /nn)nag is an 
order of magnitude lower than in the diffuse interstel- 
lar medium. In a recent study of the thermal structure 
of B68, Bergin et al. (2006) reach a similar conclusion, 
for different, physical rather than chemical reasons. Of 
course, other approaches might be adopted, for exam- 
ple, to vary the cosmic ray ionization rate, which is a 
determinant of the timescale of the ion-neutral chem- 
istry. However, test calculations have indicated that, 
unless the cosmic ray ionization rate is much larger 
than the value adopted here (( = 1 x 10~ 17 s _1 ), it is 
not as important a parameter as (ng/n^ira^; and ion- 
ization rates C >> 10~ 17 s" 1 would give rise to gas ki- 
netic temperatures which are higher than are observed. 

— The significance of the ortho:para H 2 ratio for the 
deuteration of species such as N 2 H + and NH 3 (Flower 
et al. 2006) has been emphasized. Ortho-H 2 restricts 
the degree of deuteration of H^", because the reverse 
of the reaction Hj^(HD, H 2 )H 2 D + , when it involves 
ortho-H 2 , is more energetically favourable and much 
more rapid at low temperatures than the correspond- 
ing reaction with para-H 2 . Any reduction in the abun- 
dances of the deuterated species H 2 D+, D 2 H+ and D3 



leads to a fall in the abundances of species such as 
NH 2 D, ND 2 H and ND3, which are produced in reac- 
tions with the deuterated forms of H^" . The timescale 
for the deuteration of N 2 H + , through reaction © and 
the analogous reactions with D 2 H + and D^", becomes 
comparable with the free-fall time. In the case of the 
prestellar core L1544, we have shown that the re- 
sults of our model are consistent with the observed 
high levels of deuteration of N 2 H + and NH3 when 
0.5 <J a g <J 1.0 /im. Our previous study, of freeze-out 
and coagulation during protostellar collapse (Flower et 
al. 2005), indicates that such large grain sizes cannot 
be attained through coagulation during free-fall col- 
lapse; significant coagulation of the grains must have 
occurred prior to collapse. 
— Singly deuterated ammonia, NH 2 D, was considered 
as a special case. It is special, in that it is the only 
species which has been observed in prestellar cores 
in both its ortho and para forms. As such, it pro- 
vides a test of the statistical procedures which we have 
adopted to calculate the number densities of molecules 
in specific nuclear spin states. For NH 2 D, we pre- 
dict an ortho:para ratio of 2:1, which is consistent, 
to within the error bars, with results for 4 of the 8 
prestellar cores observed by Shah & Wootten (2001). 
The deuterated forms of ammonia attain their maxi- 
mum fractional abundances for densities in the range 
5 x 10 5 ^ hh ^ 5 x 10 6 cm" 3 . The maximum fractional 
abundance shifts to higher densities with increasing de- 
gree of deuteration and with increasing grain size. In 
particular, ND3, is present with significant fractional 
abundance only in high density gas, where nn ~ 4x 10 6 
cm~ 3 . 

Appendix A: Nuclear spin statistics 

In the calculations reported above, the abundances not 
only of the chemical species but also of their individual 
nuclear spin states have been determined, where applica- 
ble. The cases of interest here involve molecules or molec- 
ular ions comprising two or more protons, such as H 2 . 
Fermi-Dirac statistics require that the nuclear wave func- 
tion should be asymmetric under exchange of identical 
protons, and this restriction associates states of a given 
nuclear spin symmetry with states of an appropriate rota- 
tional symmetry. In the case of H 2 , states with total nu- 
clear spin I — (para-H 2 ) are associated with rotational 
states with even values of the rotational quantum num- 
ber, J; states with 1=1 (ortho-H 2 ) are associated with 
odd J. Transitions between states of differing nuclear spin, 
1 = and / = 1 in this example, are induced by proton- 
exchange reactions with the ions H + , Hj, H 2 D + , D 2 H + 
and HCO + . In view of the significance of H 2 in the chem- 
istry of molecular clouds, it might be anticipated that its 
ortho:para ratio would influence the analogous ratios in 
other molecules and molecular ions. 

In order to provide a framework for the interpretation 
of the numerical results which are presented in this pa- 
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per, we establish the values of the relative populations 
of the nuclear spin states of a number of key species, in 
static equilibrium. The introduction of the distinction be- 
tween ortho and para forms of molecules and molecular 
ions leads to a large increase in the number of chemical re- 
actions which has to be considered. It is desirable but not 
feasible to consider each reaction in detail. Consequently, 
a simpler, statistical approach was adopted. 

A.l. NH%, NH3, NHl 

Let us consider the nitrogen bearing ions which are pro- 
duced in the hydrogenation sequence which commences 
with N + (ortho-H 2 , H)NH + ; the corresponding reaction 
with para-H 2 is endoergic by almost 170 K and negligible 
at T 



10 K. NH 2 forms in the reaction of NH + with H2 



NH 



H 2 



NH, 



H 



(A.l) 



for which the rate coefficient was taken to be 1.27 x 10 
cm 3 s _1 . The rate coefficients for the reactions with para- 
and ortho-H2 individually were deduced on the basis of 
the following arguments. 

— When NH + reacts with para-H2, in which the proton 
spins are "anti-parallel" , the probability that it will 
pick up the proton with a spin parallel to that of its 
own proton is the same as the probability that it will 
pick up the proton with an anti-parallel spin. Thus, we 
adopted a rate coefficient of 6.35 x 10 -10 cm 3 s~ 1 for 
the production of both ortho- and para-NH^ . 
In this case, a more rigorous argument, based on con- 
siderations of the density of states, supports the out- 
come of the simple approach. In the reaction 



NH++H 2 (p) ^NH+(p)+H 



(A.2) 



the initial and final value of the total proton spin, I, 
of the reactants and of the products) is I = 1/2. On 
the other hand, in 



NH++H 2 (p) ^NH+(o) + H 



(A.3) 



the final spin is I = 1/2 or I = 3/2. However, only 
those channels with 1=1/2 can contribute, because I 
must be conserved. Thus, although the final density of 
states, 5Z.r(2/ + 1) = 6 in this case, i.e. 3 times greater 



than in the reaction l|A.2|l . which yields para-NH 2 , 
only 2 of the 6 final states are available. Hence, the 
probability of forming ortho-NH^ is the same as that 
of forming para-NH^" , as concluded above. 
When NH + reacts with ortho-H 2 , in which the proton 
spins are "parallel" , half of the collisions will occur 
with the spin of the proton of NH + parallel to those 
of H 2 , and half with the spin anti-parallel. Hence, the 
rate coefficients for producing ortho- and para-H 2 are, 
once again, 6.35 x 10 -10 cm 3 s . 
Now consider 



and 



NH+ + H 2 (o) -> NH+(o) + H 



(A.5) 



In the first of these two reactions, I = 1/2 or I = 
3/2 for the reactants, whereas 1 = 1/2 only for the 
products. In the second reaction, I = 1/2 or I = 3/2 
for both the reactants and the products. One would 
conclude that the probability of forming ortho-NH^ 
in (|A.5(1 is 3 times that of forming para-NH^ in (|A.4|) . 
It is instructive to pursue the reasoning somewhat fur- 
ther than considerations of the density of states only. 
In the reaction (|A.5|I . when I = 1/2, ortho-NH^ may 
be considered to form by capture of N + by ortho-H 2 , 
rather than capture of a proton by NH + . If, in fact, 
N + capture is improbable for other (chemical) rea- 
sons, then the 1=1/2 channels are eliminated from 
ortho-NH^ formation. One would then conclude that 
the probability of forming ortho-NH^ is twice that 
of forming para-NH^ . Essentially, the reaction form- 
ing para-NH^ goes through an intermediate complex 
of para-NHjj~ (I = 1/2), whereas the reaction form- 
ing ortho-NH^ goes through an intermediate complex 
of ortho-NH^ (7 = 3/2). However, as in the case 
of ortho- and para-NH3, whilst the ratio of the nu- 
clear spin statistical weights is 2:1, there are twice as 
many rotational states of para as of ortho symmetry. 
In other words, the overall ratio of statistical weights is 
1:1. Thus, one concludes finally that the probability of 
forming para-NH^ , through an intermediate complex 
of para-NH^, is the same as the probability of form- 
ing ortho-NH^", through an intermediate complex of 
ortho-NH|. 

Ortho- and para-NH^ are destroyed, at the same rate, 



by para-H 2 , 

NH+(o) + H 2 (p) -> NH+(o) + H 
NH+(o) + H 2 (p) ^NH+(p)+H 
NH+(p)+H 2 (p) ^NH+(p)+H 



(A.6) 
(A.7) 
(A.8) 



NH+ +H 2 (o) -> NH+(p) +H 



(A.4) 



and hence ortho- and para-NH^ are predicted to have the 
same abundance. The rate coefficients for reactions l|A.6ll 
and (|A.7(1 are taken to be equal to one half that for reac- 
tion (|A.8(1 . and so the overall branching ratio in favour of 
para- NH+ is 3:1. Both para- and ortho-NH^ are removed 
by para-H 2 , yielding NH4" in the nuclear spin states 1 = 0, 
I = 1 or I = 2. Allowing for the 3:1 para:ortho NH^ ra- 
tio, the 1 = 0:1 = 1:1 = 2 density ratios expected 
statistically are 3:4:1, when NH^~ forms. 

According to Friedrich et al. (1981), the lowest 1 = 0, 
1 = 1 and I = 2 of NH|" are separated by only 148 mK 
in the "librational ground state". In this case, proton- 
exchange reactions with ortho-H 2 can interconvert the nu- 
clear spin states at low temperatures, but the correspond- 
ing reactions with para-H 2 are endoergic (by 170.5 K) and 
are negligible at low T. The 1 = 0:1=1:1 = 2 ratios 
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which are expected, in equilibrium, from proton exchange 
alone are, on the basis of statistical arguments, 3:4:1, i.e. 
the same as in the process of formation of NH^, consid- 
ered above. Thus, proton-exchange reactions would not 
be expected to modify the 7 = 0:1=1:1 = 2 ratios of 
NH+. 

The experiments of Friedrich et al. (1981) were per- 
formed on NH^" ions in a lattice, i.e. in the solid state, and 
their relevance to ions in the gas phase is not self-evident. 
In his book on Infrared and Raman Spectra, Herzberg 
(1945) discussed the case of CH4, which is isoelectronic 
with NH4 . On the basis of his discussion, we should con- 
clude that the nuclear spin state with the lowest energy, 
I = 2, is associated with the rotational state J = (i.e. 
the ground rotational state), I — 1 with J = 1, and 1 = 
with J = 2. Taking Ej = BJ{J+ 1) and B « 5 cm' 1 (the 
approximate value of the rotational constant of CH4; that 
of NHj" should not be much different), all from Herzberg's 
text, gives an overall separation of the nuclear spin states 
E(I = 0) - E(I = 2) « 43 K. As this separation is much 
less than the difference in the energies of the ortho and 
para ground states of H2 (170.5 K), it remains true that 
proton-exchange reactions with ortho-H2 but not para- 
H2 can interconvert the nuclear spin states of NH^ at low 
temperatures. 

A. 2. NH 3 

Dissociative recombination of NH^ with electrons yields 
both ortho-NH 3 (I = 3/2) and para-NH 3 (I = 1/2); the 
lowest state of para-NH3 lies 23.4 K above the lowest state 
of ortho-NH3. Allowing for the relative weightings of the 
population densities of the 1 = 0:1=1:1 = 2 states 
of NH4 (see Appendix IA.1(I , and the adopted values of 
the dissociative recombination coefficients, the expected 
para:ortho NH3 ratio is 3:1. The para:ortho NH3 ratio can 
be modified by proton-exchange reactions with H^, H2D + 
and D 2 H + ; but these reactions are insufficiently rapid to 
significantly affect the ratio, even in the limit of static 
equilibrium ('steady-state'), when the computed value is 
close to 3. 

NH^, the precursor of NH3, forms in reactions with 
para-~ii2- On the other hand, NH^ undergoes proton- 
exchange reactions with ort/io-r^. Consequently, we 
might expect the para:ortho NH3 ratio to depend on 
the ortho:para H2 ratio. However, as we have seen in 
Appendix IA.1I the 1 = 0:1 = 1:1 = 2 population 
density ratios of NH4 are predicted to be 3:4:1, indepen- 
dent of whether the formation or the proton-exchange re- 
actions dominate. Consequently, the para:ortho NH3 ratio 
is predicted to be independent of the ortho:para H2 ratio. 

A.3. NH 2 D 

Finally in this Section, we consider the only species 
(NH2D) whose ortho:para abundance ratio has been ob- 
served directly in prestellar cores (Shah & Wootten 2001). 
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Singly-deuterated ammonia, NH2D, is produced by the 
dissociative recombination reaction 

NH 3 D+ + e" -> NH 2 D + H (A.9) 

Its precursor, NH3D + , is formed principally in the 
deuteron-transfer reaction 

NH 3 + H 2 D+ -> NH 3 D+ + H 2 (A.10) 

and the analogous reactions with D2H+ and D3". 
Accordingly, NH3D+ is formed with the same para:ortho 
ratio as NH3, i.e. 3:1, where the lowest para state is en- 
ergetically higher than the lowest ortho state. However, 
just as in the case of NFFj", proton-transfer reactions 
with ortho-H2 modify the para:ortho ratio of NH3D+. 
Statistical considerations suggest that the rate coefficient 
for the reaction of ortho-H 2 with para-NH3D + , forming 
ortho-NH 3 D+, 

NH 3 D+(p) + H 2 (o) -» NH 3 D+(o) + H 2 (p) (A.ll) 

is the same as for the reaction with ortho-NH3D + , forming 
para-NH 3 D+, 

NH 3 D+(o) + H 2 (o) -» NH 3 D+(p) + H 2 (p) (A.12) 

Furthermore, these reactions are sufficiently rapid to de- 
termine the para:ortho NH3D + ratio, which consequently 
becomes 1:1. Dissociative recombination of NH3D+ with 
electrons favours the production of ortho-NH 2 D, relative 
to para-NH 2 D, by a factor of 2. Thus, the ortho:para 
NH2D ratio which is expected on the basis of these ar- 
guments is 2:1. 
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Fig. 4. The fractional abundances, relative to nn, of (a) 
NH3 and (b) N2H" 1 " , predicted by the free-fall model with 
an initial density tih = 10 4 cm" 3 , for grain sizes in the 
range 0.05 < a g < 1.0 fim. 5(N) = 5(0) = 0.1 and 5 = 
1.0 for all other atomic and molecular species. 
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Fig. 5. The calculated and observed values of the frac- 
tional abundances, relative to Uh, of CO, NH3 and N2H + , 
and the corresponding values of the relative abundances 
of the deuterated forms of NH3 and N2H+. The calcu- 
lated values derive from the free-fall model with an initial 
density tih = 10 4 cm~ 3 and grain radii a g = 0.5 /im and 
a g = 1.0 /xm. 5(N) = 5(0) = 0.1 and 5 = 1.0 for all 
other atomic and molecular species. The observed ratios 
are denoted by the crosses, with error bars. 



Fig. 6. (a, c, e) The abundance ratios of the deuterated 
forms of N2H + and NH3, together with (10 4 times) the or- 
tho:para H2 ratio, and (b, d, f) the fractional abundances 
of the deuterated forms of ammonia, in the gas phase, as 
functions of rin; (a) and (b) correspond to a g = 0.2 /xm, 
(c) and (d) to a g — 0.5 /xm, (e) and (f) to a g = 1.0 /J,m. 
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Fig. 7. (a) The fractional abundances of and its 
deuterated isotopes, and (b) the relative abundances of 
their ortho and para forms (in the sense of excited state 
divided by ground state abundance), for a g = 0.5 /una. The 
ortho:para H2 ratio is also plotted, for reference. 



